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Entropy of a Quantum Oscillator coupled to a Heat Bath and 
implications for Quantum Thermodynamics 

G. W. Ford a , R. F. 0'Connell b 
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The free energy of a quantum oscillator in an arbitrary heat bath at temperature T is given by a "remarkable 
formula" which involves only a single integral. This leads to a corresponding simple result for the entropy. The 
low temperature limit is examined in detail and we obtain explicit results both for the case of an Ohmic heat 
bath and a radiation heat bath. More general heat bath models are also examined. This enables us to determine 
the entropy at zero temperature in order to check the third law of thermodynamics in the quantum regime. 



1. Introduction 

The widespread interest in recent years in (a) 
mescoscopic systems |1I2I3I4I5I6I7| and (b) funda- 
mental quantum physics and quantum computa- 
tion |8l9ll0lilll2ll3ll4ll5llfi| has highlighted the 
critical role which dissipative environments play 
in such studies. This has led to a critical ex- 
amination of many results that were derived for 
macroscopic systems. In particular, there has 
been considerable interest in the area of quan- 
tum and mesoscopic thermodynamics, the sub- 
ject of this conference. In particular, in some in- 
stances questions have been raised about the va- 
lidity of the fundamental laws of thermodynam- 
ics. Whereas many interesting new facets of old 
results have emerged, it is important to exercise 
caution before questioning the validity of funda- 
mental laws (especially the three laws of thermo- 
dynamics), since many subtle issues arise. 

Here, we examine the third law of thermody- 
namics in the quantum regime by calculating the 
entropy S for a quantum oscillator in an arbitrary 
heat bath at temperature T and checking to see if 
it vanishes as T — > 0, in conformity with Nernst's 

law [nj. 

The first question which arises is how to calcu- 
late S. The von Neumann result for S involves 
p log p, where p is the density matrix for the whole 
system, a non-trivial quantity to calculate. As 



a result, one's first thought is perhaps to make 
use of the Wigner distribution function W corre- 
sponding to p. This has led some authors to sim- 
ply replace p with W in the von Neumann result 
but, unfortunately this is not correct because the 
W corresponding to log p is not the log of the W 
corresponding to p and, as a consequence, one is 
led to misleading conclusions. A way out of this 
impasse is to use the method which we introduced 
in 1985, in collaboration with J.T. Lewis ^Hjj in 
order to calculate the free energy F. Then, us- 
ing a familiar thermodynamic relation, the result 
for S readily follows. In addition, the result for 
the total energy U follows in a similar manner. 
In Section 2, we review this method and write 
down the results for F, S, and U in the case of 
a quantum oscillator in an arbitrary heat bath at 
temperature T. All of these results involve just 
a single integral. In section 3, we evaluate the 
relevant integral for the case of an Ohmic heat 
bath and for an arbitrary T. This enables us to 
show that S — f as T — > 0, in conformity with 
Nernst's law. A similar result is obtained in Sec- 
tion 4 for the case of a blackbody radiation heat 
bath and in Section 5 for an arbitrary heat bath. 
Our conclusions are given in Section 6. 
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2. Fundamentals 

A very general model for the motion of a quan- 
tum particle in an arbitrary heat bath is the 
independent-oscillator model |14I18| , which is de- 
scribed by the Hamiltonian 



H 




x y 



(1) 



Here m is the mass of the quantum particle while 
Tfij and ujj refer to the mass and frequency of 
heat-bath oscillator j. In addition, x and p are 
the coordinate and momentum operators for the 
quantum particle and qj and pj are the corre- 
sponding quantities for the heat-bath oscillators. 
The infinity of choices for the m,j and u>j give 
this model its great generality. Moreover, as em- 
phasized in ^3], the most general coupling of a 
quantum particle to a linear passive heat bath is 
equivalent with an independent oscillator model, 
which is described by the Hamiltonian given in 

CD- 

Use of the Heisenberg equations of motion leads 
to the quantum Langevin equation 



mx + I dt'/j,(t-t')x(t')+V'(x) 

J — OO 

= F{t) (2) 

where V'(x) = dV(x)/dx is the negative of the 
time-independent external force and /i(t) is the 
so-called memory function. F(t) is the random 
(fluctuation or noise) operator force with mean 
(F(t)) = 0. The quantities fj,(t) and F(t) describe 
the properties of the heat bath. 

In the particular case of an oscillator potential 



V(x) = —Kx 2 = —muJnX 2 . 
w 2 2 



(3) 



Substituting ® into J5J) enables us to obtain the 
explicit solution 



x(u>) = a(u})F(ui), 



(4) 



where the superposed tilde is used to denote the 
Fourier transform. Thus, x(u>) is the Fourier 



transform of the operator x(t): 

POO 

x(w) = / dtx(t)e luJt , 



(5) 



and similarly for F(u>). Here ct(z) is the familiar 
response function (generalized susceptibility) 

a i z ) — 5 . —, i —■ (6) 

v ; -mz 2 - izfi(z) + K w 

and ji{z) is the Fourier transform of the memory 
function: 



fi(z) = / dt f i(t)e l 



(7) 



We have now all the tools at our disposal nec- 
essary to obtain thermodynamic qualities for the 
heat bath. Our main task will be the calculation 
of the free energy F, which is a thermodynamic 
potential from which other thermodynamic func- 
tions can be obtained by differentiation. The en- 
tropy is the only one of interest here and is given 
by the relation 

— (S)/ 

The system of an oscillator coupled to a heat 
bath in thermal equilibrium at temperature T has 
a well-defined free energy. The free energy as- 
cribed to the oscillator, F(T), is given by the free 
energy of the system minus the free energy of the 
heat bath in the absence of the oscillator. This 
calculation was carried out by two different meth- 
ods |18ll9j leading to the "remarkable formula" 

1 f x , ( d loga(w + i0+ 



F{T) 



djjjf{uj, T)Im 



io \ dm 

where f(ui,T) is the free energy of a single oscil- 
lator of frequency u, given by 

f(uj,T) = fcTlog[l -exp(-tiu/kTj\. (10) 

where the zero-point contribution (fvjj/2) has 
been omitted. Thus, all that remains is to specify 
ji(z) which characterizes the heat bath. In the re- 
maining sections, we consider various heat bath 
models. In the low temperature case (kT << 
hu>o), explicit results may be obtained by not- 
ing that f(u>,T) vanishes exponentially for large 
(Jko/kT). Hence, the integrand in Q is confined 
to small (Hui/kT) so that the factor multiplying 
f(to, T) can be expanded in powers of ui. 



(9) 
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3. Ohmic heat bath 

This is an oft-studied model for which 

#H=m7, (11) 

where 7 is a constant. Thus 

a(uo) = \m (uo 2 — uo 2 — iwy)] 1 (12) 

is the familiar phenomenological Drude-Lorentz 
model result. Hence, using lfT2T) and lfTU|) in (J5J, 

we obtain 

kT 



F(T) = / °° duo log [1 - exp (-Hw/kT)] 



Hence, in the low temperature case, 



(13) 



Im 



f dlogq(^) 1 



7(" 2 +" 2 ) 
( W 2 -^ 2 )"+ 7 2 ^ 2 



(14) 



Substituting into l|13|) and changing the variable 
of integration to 



y 



dylog(l-e^) 



(15) 



(16) 



we obtain 

^ - ^ r 

The following integral is relevant (now and later) 
dyy z log(l - e-") = -r(« + l)£(a: + 2), (17) 

where £ is the Riemann zeta-fvmction, 
00 1 

= £ -■ ( 18 ) 

71=1 

If z is an even integer £(2) is related to the 
Bernoulli numbers, C(2) = tt 2 5i = tt 2 /6, C(4) = 
7r 4 i?2/3 = 7r 4 /90, etc. But in Section 5 other val- 
ues appear. Thus, in particular, 



dy dog (1 - e 



dyy 2 log(l-e y ) 



7T 

45' 



(19) 
(20) 



Hence 

F(T) = 
Also 



7T^ / fcT 

6 V^o 



<9F 



<9T 



3 ' huoi 



(21) 



(22) 



We note that S(T) as T -> 0, in confor- 
mity with the third law of thermodynamics. 

4. Blackbody radiation heat bath 

In this case, we obtained [23j 

a(oj) = (1 - iuoT e )a D (uo), (23) 

where aofw) is the Drude result with 7 replaced 
by 



7e = 

and 



2e 2 
3mc 3 



6 x l(T 24 s. 



(24) 



(25) 



Thus, proceeding as in the Ohmic case and once 
again letting uo — > in order to calculate the re- 
sults for small T, we obtain HT\ 



Im 



d\oga(w) 
duo 



7e (w 2 + uol) 



le 



{U0l ~ + lW 1 + 

(26) 



37e 2 



It follows that 

le(kTf. 



F(T) 



nhuo 2 



kT 

huOn 



tP_ ( kT 
15 e \TujOq 

from which it follows that 

4^ 2 k 4 T 3 
b ~ 15 7e (^ ) 4 ' 



dyy 2 log(l - e y ) 



(27) 



(28) 



Once again, we see that S — » as T — > 0, in 
agreement with Nernst's law. In this case, S ^ 
faster than in the Ohmic case, as a result of the 
fact that the right-side of (|25|l has a factor uo 2 
whereas the corresponding result on the right-side 
of (|14f> is independent of uo. 
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5. Arbitrary heat bath 

From © and 10, we obtain 
Im a(uj) — uj\a(uj)\ 2 Re fi(u>), 

Re a{uj) — {— m (oj 2 — lu 2 ) 

—oj Im jl(oj)} \a(oj)\ 2 . 



(29) 



(30) 



Using these results in (0 leads to 

F(T) = i / °° duffa T)\a\ 2 {m (uj 2 + uj 2 ) 

Re jj,(uj) — uj 2 Re jl{uj)-^ Im £l(u)) 
+ lu (—moj 2 + muj 2 + ujlm fi) ^;Re /i(w)} -(31) 

Now we make use of the fact that fl(z) must be a 
positive real function |14j and hence the boundary 
value of fl(z) on the real axis has everywhere a 
positive real part i.e. 



Re [ji (u + i0 + )] > 0, - oo < uj < 



(32) 



Thus, in the neighborhood of the origin, a very 
general class of models are incorporated by writ- 
ing 



fi(z) = mb 1 - a {-iz) a , -!<«<!, 



(33) 



where b is a constant with the dimensions of fre- 
quency. It is easy to verify that this is a positive 
real function if and only if a is within the indi- 
cated range [we choose the branch — it < 8 < it 
where 9 is arg(z)]. Hence 



Re jiluj) - cu a cos {a (6 - vr/2)} 
Im p,(u) - uj a sin {a (9 - tt/2)} . 



(34) 



As it turns out, the case a = — 1 requires special 
treatment so we will consider this first, obtaining 
fl = imb 2 /z, so that 



Ot{<jj) 



1 



-mio 2 + m(b 2 + Wq) 



(35) 



Thus, this corresponds to a shift in the force con- 
stant K with no dissipation. In the absence of 
dissipation 



.d\oga{uj + i0+) ' 

Im < ; > = 7T [d(uj — U>0) 



+S(u + Wo)] 



(36) 



F(T) = f((J ,T) ^ -kTe- hu)/kT , 



(37) 



from which it is clear that S(T) — > as T — > 0. 
We now return to the more general case — 1 < a < 
1 and again consider the low-temperature limit. 
Thus, as u> — > 0, the { } term in (|3l]l) reduces to 



{ } — > too;q(1 + a)Re fl, 

and, in addition, |a(ti;)| 2 
that 



F(T) 



where 



kT 



nJ 2 ) 



(38) 
It follows 



}/H,(39) 



= r o °° dujuj a log [1 - exp (-huj/kT)] 
= m a+1 f Q °°dyy a \o g (l-e-y). 

It follows from (jSHJ, $M and © that 
xT a+1 T(a + l)£(a + 2). 



(40) 



(41) 

= 7, 



As a check, we note that for a = and b 
this result reduces to the Ohmic result given in 
(J23- Hence, since (a + 1) can never be negative, 
we conclude that 



S(T) - (kT) a+1 - 
6. Conclusions 



as T — > 0. 



(42) 



For the case of a quantum oscillator coupled to 
an arbitrary heat bath, we found that Nernst's 
third law of thermodynamics is still valid. 
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